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Experimental generation of single photons via active multiplexing 
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An on-demand single-photon source is a fundamental building block in quantum science and technology. We 
experimentally demonstrate the proof of concept for a scheme to generate on-demand single photons via actively 
multiplexing several heralded photons probabilistically produced from pulsed spontaneous parametric down- 
conversions (SPDCs). By utilizing a four-photon-pair source, an active feed-forward technique, and an ultrafast 
single-photon router, we show a fourfold enhancement of the output photon rate. Simultaneously, we maintain 
the quality of the output single-photon states, confirmed by correlation measurements. We also experimentally 
verify, via Hong-Ou-Mandel interference, that the router does not affect the indistinguishability of the single 
photons. Furthermore, we give numerical simulations, which indicate that photons based on multiplexing of 
four SPDC sources can outperform the heralding based on highly advanced photon-number-resolving detectors. 
Our results show a route for on-demand single-photon generation and the practical realization of scalable linear 
optical quantum information processing. 



Numerous applications for on-demand single-photon 
sources have been proposed in the field of quantum-enhanced 
science and technology [1-6|. The commonly studied meth- 
ods to generate on-demand single photons are based on emis- 
sions from molecules 00, atoms J9|, color centers in dia- 
mond 1 10 1, quantum dots lfTTl[T2l . and donor impurities |[T3l . 
However, each of these methods has certain challenges to 
overcome. For atomic systems, the repetition rates, the collec- 
tion efficiency, and the complexity of the experimental setup 
are the main obstacles. For color centers, quantum dots and 
donor impurities, it is difficult to achieve a time-bandwidth 
limited emission spectrum lfl4l and indistinguishable output 
photons from different sources lfT31 . 

An alternative approach is the so-called heralded single- 
photon source (HSPS) based on the process of spontaneous 
parametric down-conversion (SPDC), where a pair of photons 
is created from a nonlinear crystal pumped by a laser. The 
detection of one trigger photon indicates, or heralds, the pres- 
ence of its twin. Because of the phase-matching condition and 
energy conservation, a HSPS from SPDC has some appealing 
features. The linear momentum (direction), polarization, and 
wavelength of that heralded photon are well defined via the 
measurement of its twin. The scheme faces two main chal- 
lenges: (1) The conversion process is inherently random and 
hence there is no prior knowledge of when the heralding event 
will occur. (2) Because of the nature of the SPDC process, in 
addition to the probability of generating one pair of photons, 
Pi, there is also a finite probability of generating more than 
one pair, P>i, via higher order emissions, which decreases 
the quality of the HSPS. This latter generation probability in- 
creases nonlinearly with the interaction strength of the pump 
laser and the nonlinear crystal. 



To surmount the random production of photon pairs from 
SPDC, a pulsed laser can be used to pump the nonlinear crys- 
tals, and hence the photon pairs can be created only at certain 
times. When the pulse duration of the pump laser is much 
shorter than the coincidence measurement time, the probabil- 
ity of generating n pairs of photons per pulse can be repre- 
sented as a Poisson distribution: 



P, 



N n e 



n a -N 



forn = 0,1,2,... 



(1) 



where the mean number of pairs, N, depends on the pump 
power and parameters of the crystals. Note that it has been 
shown that to use Bose (thermal) distribution, the analysis pre- 
sented below would yield essentially the same results [ 17). An 
ideal HSPS would have exactly one photon generated from 
each pulse, which requires N = 1. However, for N = 1, 
the probability to observe more than one photon is 0.26 due to 
the Poisson distribution. This significant higher order emis- 
sions probability, P>i, drastically reduces the quality of the 
HSPS. For this reason, pulsed SPDC is normally operated at 
low power (thus P yl <1) and obeys a trade-off between the 
count rate and quality of the HSPS. This is the reason why 
a single SPDC source is fundamentally limited. This limi- 
tation can be overcome by using spatial or temporal multi- 
plexing of several sources [16-20|. Although this idea has 
attracted significant attention, as highlighted in recent review 
articles ll2Tll22l . it has not yet been experimentally realized. 
Here we follow and extend the proposal in Ref. [16]. We re- 
port an experimental realization of a 4-SPDC single-photon 
multiplexing system and demonstrate a proof-of -principle en- 
hancement over a 1-SPDC single-photon source. 



Theoretical analysis 
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In the spatial multiplexing scheme, m SPDC sources are 
pumped by a single pulsed laser whose power is split equally 



for each SPDC source by a series of beam splitters (BSs). 
These SPDC sources are then coupled by fast photon routers 
and directed to a single output. Each individual SPDC source 
has the probability to generate one pair, Pi . The pump power 
input is chosen low enough that the generation of more than 
one pair, P>i, is much smaller than Pi. With sufficiently 
large m and feed-forward operation of the fast photon router, 
the probability of obtaining one pair emission in this array, 
Qi, can approach unity and hence be on demand. The de- 
tection signal of the heralding (trigger) photon of an individ- 
ual source is used to identify which source has produced a 
pair of photons and to control the photon routers to direct 
the successfully created photon to the single output of this 
771-SPDC array. Shapiro et al. proposed a modular configu- 
ration by using 2x1 routers [16|, where two inputs can be 
routed to a single output. As shown in Fig.[TTa), in a 1-SPDC 
HSPS, one nonlinear crystal (NLC) cut for collinear type-II 
phase matching is pumped by laser pulses, where a pair of 
orthogonally polarized photons is generated. Their quantum 
state is a product state in the polarization degree of freedom: 
I$)i2 = |V), \H) 2 . |V)j and \H) 2 denote the vertical and 
horizontal polarization states of photons 1 and 2 respectively. 
Then they are separated by a polarizing beam splitter (PBS). 
The detection of the vertically polarized photon heralds the 
existence of the horizontally polarized photon, which is the 
output of this 1-SPDC source. Two such 1-SPDC sources can 
be coupled with a photon router and integrated into a 2-SPDC 
module, as also shown in Fig. ITTa). It is straightforward to 
increase the size of the network. Two 2-SPDC modules can 
build a 4-SPDC module and so on [Fig. [Tib)]. The advantages 
of this multiplexing are as follows: First, one can enhance the 
one-pair-generation probability by a gain factor of (see Ap- 
pendix A) 
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Second, the signal-to-noise ratio of an m-SPDC array is the 
same as that of a 1-SPDC source: J^- — -J^-. Q^i is the 
probability of obtaining more than one pair in this array. In 
Fig.fljc), the gain G of a multiplexed system with m sources 
(blue) is shown and compared with a single 1-SPDC source 
(red) at N = 0.1. For P close to 1 (N close to 0), it first 
increases linearly with m. The gain saturates at the value 
limm^ooG = jjp = 1 _ 1 _ JV for an infinite number of 
sources. The black dashed line in Fig. [TJc) shows this satu- 
ration limit at N = 0.1. The dependence of the gain on the 
mean photon pair number per SPDC source N and the num- 
bers of the sources m is shown in Fig. [JJd). Intuitively, in 
order to obtain high-quality and high-rate single photons, one 
should operate in a regime, in which N is low for each SPDC, 
in order to reduce the higher order emission per source, and 
77i is large, so that the chance of one source firing in this array 
is high. This is the reason why the optimal operating com- 
bination of N and m in terms of maximizing G lies in the 
upper-left corner in Fig. [JJd), given that m is limited to 10. 
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FIG. 1 : (Color online) Scheme and calculations of generating on- 
demand single photons via multiplexing with active feed forward 
from spontaneous parametric down conversion (SPDC). (a) A 2- 
SPDC module consists of two 1-SPDC sources coupled with a pho- 
ton router. Each 1-SPDC source produces a pair of photons proba- 
bilistically from a nonlinear crystal (NLC) pumped by laser pulses 
whose power is split by a beam splitter (BS). The signal of each 
avalanche photodiode (APD) detector, indicated by the black arrows, 
is feed-forwarded and controls the photon router to direct one of the 
successfully created photons to a single output. This output single 
photon along with the output feed-forward signal can then be used 
in the next level. Suitable individual optical delays (not shown) are 
incorporated in order to erase the output photons' temporal distin- 
guishability from different sources, (b) Similarly, a 4-SPDC mul- 
tiplexing system can be realized by coupling two 2-SPDC modules 
with another photon router. An m-SPDC multiplexing system re- 
quires m SPDC sources and (m — 1) routers, (c) The calculated 
gain, G, of using a multiplexed system is plotted as a function of the 
number of SPDC sources, m, where only integer values of m is rel- 
evant (blue solid line). The trivial gain of using a 1-SPDC source is 
shown for comparison (red solid line). The mean photon pair number 
per SPDC source is N = 0.1, indicated with the vertical dashed line 
in (d), and the horizontal black dashed line represents the saturation 
limit lim„ woo G = — ^__ ~ 10.5. In addition, we show the gain 
when using multiplexing with a 0.95 transmission per path router and 
0.97 of the polarization router (blue dashed line) in order to illustrate 
loss effects. For details, see the main text and Appendix A. d, Calcu- 
lation result of the gain, G as a function of both the numbers of SPDC 
sources m and the mean photon pair number per SPDC source N. 



Experiment 

Initially, we created photon pairs from a /3-barium-borate 
crystal (BBO) via noncollinear type-II phase matching [23|. 
As shown in Fig. 2a, photons 1 and 2 were generated from 
BBOl and their quantum state was | <E> > I2 . = T7s{\H)\ \V) 2 — 
|V")j \H} 2 ). In comparison with the product state gener- 
ated from the collinear phase matching (\V) l \H) 2 ), this non- 
collinear scheme is equivalent to having two 1-SPDC sources 
(\H) l \V) 2 and \V) ] \H) 2 ), each subject to half of the pump 
power. Additionally, the non-collinear type-II phase matching 
overlaps the spatial modes of these two sources as schemati- 
cally shown by the dashed PBS in Fig. 2(c). By using a po- 
larization switch and a polarizer, these two 1-SPDC sources 
form a 2-SPDC module. In comparison with the scheme pro- 
posed in Ref. [16], our noncollinear scheme has two advan- 
tages: The amounts of crystals and beam splitters needed for 
building an m-SPDC array are reduced from to to ^ and there 
is no need to spatially overlap the photons generated from two 
1-SPDC sources. An identical 2-SPDC module subsequently 
generated photons 3 and 4 from BB02. See Appendix B for 
details. 

These two 2-SPDC modules form a 4-SPDC multiplexed 
system with an ultrafast photon router, which satisfies two 
main criteria: (1) short response time, which is necessary for 
the feed-forward operations in a large system, and (2) high 
routing visibility, which is needed for a high quality HSPS. As 
shown in Fig. 2(b), we employed a Mach-Zehnder interferom- 
eter (MZI) to realize the ultrafast photon router. It consisted 
of two 50:50 BSs, mirrors, and most importantly two electro- 
optic modulators (EOMs) with one in each arm of the MZI, 
where two rubidium titanate phosphate (RTP) crystals were 
used as the electro-optic material. This design of the router 
has several advantages: (1) It can route incoming photons 
independent of their polarizations, which is crucial for cer- 
tain tasks. In the schemes of polarization-entanglement-based 
quantum cryptography ll24l and quantum computation with 
polarization encoded photonic qubits |6|, any polarization- 
dependent routing will destroy the polarization entanglement 
and hence defeat the original goal of achieving communica- 
tion security and computation speed-up provided by the en- 
tanglement. For instance, in polarization-entanglement-based 
quantum cryptography described in Ref. 11241 . the source is 
located at Alice's laboratory, and out of each pair of entan- 
gled photons, one is sent to Alice and the other to Bob. In 
a multiplexed scenario, many sources produce polarization- 
entangled photon pairs. The trigger photons are sent to and 
measured by Alice. Depending on the outcome of her po- 
larization measurement, she sends an electronic signal to the 
routers for multiplexing the photons which will be sent to Bob. 
Multiplexed systems with our router can enhance the key rate 
while keeping the quantum bit error rate (due to the multi- 
ple pair generations) constant. Whatever critical key rate is 
needed at the receiver, to-SPDC multiplexing extends the pos- 
sible communication distance between the two parties by the 
additional distance ln(m)Lo compared to 1-SPDC, where L$ 
is the 1/e decay distance in the transmission channel (valid for 
both free space and fiber). In a 1550 nm fiber-based cryptogra- 
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FIG. 2: (Color online) Experimental setup, (a) Two pairs of the po- 
larization entangled photons (1 and 2, 3 and 4) were generated from 
two /3-barium borate crystals (BBOl and BB02) via SPDC, where 
half-wave plates (A/2) and compensating BBO crystals were used to 
counter walkoff effects in the down-conversion crystal, and single- 
mode fibers and interference filters (IFs) were used to clean their 
spatial and spectral modes, respectively. Photons 1 and 4 entered 
into spatial modes a and d respectively. Photons 2 and 3 entered 
into spatial modes b and c and were delayed in single-mode fibers 
(~ 97m) with respect to photons 1 and 4. This allowed the logic and 
electronics to have enough time to control the EOMs in the router. 
Due to the sequential configuration of the sources (see Appendix A), 
appropriate individual fiber delays are implemented to eliminate the 
output photons' temporal distinguishability from different sources. 
The layout of the ultrafast photon router is shown in (b). See text and 
Appendix Cfor details and notations, (c) The experimental scheme 
of a 4-SPDC multiplexed system with active feed-forward. Sources 
I and II, III and IV generated photon pairs from BBO 1 (BBO 2). 
Sources I and II (III and IV) were inherently overlapped via the non- 
collinear phase matching, as explained in the text. The polarization 
states of photons 1 and 4 were measured with polarizing beam split- 
ters (PBS) and single-photon counting modules (SPCMs). The de- 
tection signals of photons 1 and 4 were fed into the logic and then 
used to control the ultrafast photon router. Depending on whether 
photon 1 or photon 4 was detected, the ultrafast photon router would 
be switched on or off. Thus, all the photons exited from mode b" and 
no photons exited from c". The quality of the output single photons 
was quantified via measurements of the correlation function, <? 2 (0). 



phy system, a 4-SPDC multiplexing system could increase the 
available distance by ln(4)21.7km w 30.1 km. Note however, 
whether one needs entangled photons depends on what kind of 
quantum-information processing task the photon will be used 
for later on. One does not need entanglement for the oper- 
ation of the multiplexing procedure itself or for using mul- 
tiplexing in the Bennett and Brassard quantum cryptography 
protocol 1 25 1, as the latter is not based on entanglement. Then 
one can also use a mixture of polarization-correlated product 
states as the input of the multiplexing scheme. (2) A simi- 
lar design has been widely used in the photonic industry to 
implement intensity modulators |26l . Therefore, the scheme 
demonstrated here has the potential to be implemented with 



the mature technology of integrated optics. 

The optical axes of both RTP crystals in the EOMs were 
oriented along 45° and the voltages applied to them were al- 
ways of the same amplitude but with opposite polarities. If 
the opposite half-wave voltages (HWVs) were applied to the 
EOMs, a 7r phase shift would be introduced to the MZI. We 
locked the phase of the MZI at when the EOMs were off. 
This phase is defined via the condition that all the photons 
that enter from input b exit into c", and c into b" (Fig. 2b). 
Therefore, when the EOMs were on and were applied with 
opposite HWVs, the photons that entered from input b exited 
into b", and c into c". 

The rising and falling times of our router for a 7r-phase 
modulation are each about 5.6 ns. (This could be reduced to 
tens of picoseconds with state-of-art integrated photonics de- 
vices [26 1.) The recharger time of the high-voltage driver of 
our EOMs is about 50 ns. The traveling time in the used coax- 
ial cables is also a few nanoseconds. Together, this in prin- 
ciple allows a repetition rate of our router of approximately 
15 MHz. (Note that the dead time of trigger APDs is about 
150 ns, which will limit the operating frequency of small- 
scale multiplexing systems.) The visibility of 7r-phase modu- 
lation (routing visibility), V", was above 95%. It is defined 

r„-i°„ 
as V£ = X +J b for spatial mode b as the input, where Ig„ 

and IP,, are the output intensities of spatial mode b" of phase 
and 7r, respectively. Similar results were obtained for spa- 
tial mode c as the input. Therefore, this photon router meets 
the two requirements: short response time (compared to the 
recharging time of our EOMs) and high routing visibility. See 
Appendix C for further details of this ultrafast photon router. 

The 4-SPDC multiplexing system [Fig. 2(c)] was operated 
as follows. We detected photons 1 and 4 in the modes of 1H, 
IV, 4H and 4V, and sent the electrical signals from the detec- 
tors through a custom-built programmable logic to the router. 
When the detector on mode 1H or IV fired, the EOMs in the 
MZI were switched on with HWV and thus the phase of the 
MZI was switched to ir. Therefore, photon 2 was routed to 
spatial mode b". When the detector on mode 4H or 4V fired, 
the EOMs in the MZI were switched off, and the phase of the 
MZI remained at 0. Therefore, photon 3 was routed to spa- 
tial mode b". In the last step, an additional electro-optical 
modulator (EOM3) placed in spatial mode b" acted as a po- 
larization router and performed an active polarization rotation 
on the output photon, before it was finally cleaned by using a 
PBS. However, due to limited resources, we bypassed EOM3 
and performed this operation with post-selection. Note that 
introducing EOM3 would only reduce the routing visibility 
by at most 1 % (with a polarization switching visibility of at 
least 99%) EH. 

In order that our ra-SPDC source is useful in photonic 
quantum-information processing, it is most important that it 
performs well in Hong-Ou-Mandel (HOM) two-photon in- 
terference [28|, which is at the heart of many protocols and 
has been demonstrated with different systems and schemes. 
For instance, interference between consecutive photons emit- 
ted from the same quantum dots [12], between photons from 
different donor impurities [ 1 3 1 and molecules [8|, as well as 



between two independent SPDC sources 11291 1301 have been 
shown. At present, in terms of the count rate and qual- 
ity of the independent single photons, femtosecond-pulsed 
SPDC is still better than other schemes and has been success- 
fully used in multiphoton experiments in the past. A single- 
photon source based on our scheme is suitable for high quality 
HOM two-photon interference and enables multiphoton ex- 
periments [31]. The present multiplexing scheme includes the 
photon router as the only additional complexity which was 
not involved in earlier experiments. Since the path length 
difference of the MZI can be controlled very accurately, the 
router will not introduce any temporal distinguishability, even 
though it is based on first-order interference which is wave- 
length sensitive (nm range). In contrast, most of the pho- 
tonic quantum-computation experiments (C-phase gate ll32l - 
[34), entanglement swapping fl4j [5j, etc.) rely on second-order 
interference, which is coherence length sensitive (depends on 
the context, but at least more than 10 4 nm). In order to empiri- 
cally exclude possible degradations of the indistinguishability 
due to the router, we have performed HOM two-photon inter- 
ference experiments with the router. 

Experimentally, we use a pair of photons generated from 
one source and send it to the router. We vary the path- 
length difference between these two photons with a motorized 
translation stage mounted on one of the fiber coupling stages 
and measure the twofold coincidences between two detectors 
placed directly behind two outputs of the router [b" and c" 
in Fig. ETb)]. The phase of the MZI and hence the reflec- 
tivity (and transmissivity) of the router is varied by applying 
different voltages on the EOMs. In case of a ir/2 phase, the 
router becomes a balanced BS and the distinguishability of 
two input photons' spatial modes is erased. In consequence, 
the minimum of the coincidence counts occurs for the optimal 
temporal overlap (with the help of suitable individual fiber de- 
lays) of the two photons, and HOM two-photon interference 
with a visibility of 88.7% ± 3.8% can be observed (red circles 
in Fig. [3]). In case of a phase, the whole router represents 
a highly transmissive BS (routing visibility above 95%) and 
the two photons remain distinguishable in their spatial modes. 
Correspondingly, one obtains path-length-difference insensi- 
tive coincidence counts (black square in Fig. RJ. This is in 
agreement with complementarity, where in principle no HOM 
interference can be observed. In addition, we have compared 
the results when using the router with those when using a nor- 
mal free-space BS: quantitatively identical results have been 
obtained (within the error bars). These HOM experiments to- 
gether with the earlier works |29, 31] demonstrate the suit- 
ability of the multiplexed photons for scalable quantum infor- 
mation processing. 

We quantified the quality of the multiplexed single pho- 
tons with the second-order correlation function at zero time 
delay, g 2 (0). The smaller g 2 (0) is, the higher the quality of 
the single photons becomes. The correlation function can be 
measured in good approximation with a 50:50 BS and two 
more detectors 1351 . as shown in Fig. 2c. It is defined as 
5 2 (0) = N tTR N t /(N tT N tR ), where N tTR , N tT , N tR , and 
N t are the coincidence counts among trigger, transmission, 
and reflection of the beam splitter, the coincidence counts be- 
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FIG. 3: (Color online) Experimental HOM two-photon interference 
performed with the photon router. The two-fold coincidences are 
measured between two detectors placed directly behind two outputs 
of the photon router (b" and c" in Fig. pb) and are plotted over the 
relative delay of the interfering photons [realized by changing the 
length of path c in Fig. 2(b)]. By adjusting the voltage applied on 
the EOMs, we set the phase of the MZI and hence varied the re- 
flectivity (and transmissivity) of the router. When the phase of the 
MZI is 7r/2, HOM two-photon interference shows up (red circles). 
This is because the distinguishability of the two input photons' spa- 
tial modes has been erased by the router. Consequently, the mini- 
mum coincidence counts occur for the temporal overlap of the two 
photons, i.e. zero time delay. The red solid line is a Gaussian fit to 
the data with a visibility of 88.7% ± 3.8%. When the phase is 0, path 
length difference insensitive coincidence counts are obtained (black 
squares) and no HOM two-photon interference shows up, because 
the two photons remain distinguishable in their spatial modes. The 
black solid line is the average of the coincidence counts for phase 
setting. Background has not been subtracted. Error bars indicate ±1 
standard deviation. 



tween trigger and transmission, the coincidence counts be- 
tween trigger and reflection, and the single counts of the trig- 
ger, respectively. Figure H shows g 2 (0) versus the counts per 
350 seconds with 1-SPDC, 2-SPDC and 4-SPDC multiplexing 
systems. 

We performed the experiments with various levels of pump 
power (100, 250, 500, and 1000 mW) to demonstrate the her- 
alded single-photon quality and generation rate trade-off of 
the SPDC sources. With a 1-SPDC source (black circles in 
Fig.HJi, the counts increased as we increased the pump power. 
The quality of the output single photon decreased as g 2 (0) 
was increased because of higher order emissions. At a con- 
stant mean photon pair number N per source, with single- 
photon multiplexing, 2-SPDC (red triangles in Fig. HI and 4- 
SPDC multiplexed (blue squares in Fig. HI systems enhanced 
the counts by factors of 2 and 4, respectively, while g 2 (0) 
remained the same as the 1-SPDC. On the other hand, one 
can understand Fig.HJalso from the following perspective: By 
accordingly decreasing TV per source, one can employ multi- 
plexed systems while keeping the counts constant and increase 
the quality of the output single photon via decreasing g 2 (0). 
In the regime of small N, the counts scale proportionally to m 
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FIG. 4: (Color online) Experimental results. The correlation function 
at zero delay g 2 (0) is plotted versus the counts in 350 s for 1-SPDC 
(black circles), 2-SPDC (red triangles) and 4-SPDC (blue squares) 
multiplexed systems, where all of them contain the router circuitry. 
The g 2 (0) function saturates at the value 1 for large count rates, i.e. 
for large pump powers. The g 2 (0) data shown are in linear regime of 
low count rates and hence fitted with lines through the origin. Using 
multiplexed systems, one can increase the count rate while keeping 
constant the signal-to-noise ratio and hence the <? 2 (0) function. In 
other words, increasing the number of SPDC sources enables one 
to increase the count rate with the same signal-to-noise ratio or (if 
one turns down the pump power) improves the quality of the single- 
photon output while maintaining the count rate of the output. Error 
bars indicate ±1 standard deviation. 



at fixed <? 2 (0), and g 2 (0) scales inversely proportional to m at 
fixed counts. Because the coupling efficiencies of the sources 
are different, the count rates and correlation functions are also 
different. The values of the counts and the correlation func- 
tion shown in Fig. HJ were averaged over four SPDC sources 
(see Appendix B). 



Discussion 

The transmission throughput of the single-photon router is 
an important parameter for the viability of a multiplexed sys- 
tem, and next we provide an analysis of losses in our sys- 
tem. In our experiment, the output photon's transmission of 
the 4-SPDC multiplexing system is T « 0.1, where the main 
loss is due to single-mode-fiber to single-mode-fiber coupling 
in the router. To have an advantage of 4-SPDC multiplex- 
ing over a 1-SPDC without router circuitry would necessitate 
T > 1/4 or, equivalently, T > 1/2 if the routers are all 
the same and each have a transmission of Tq. The gain due 
to multiplexing is G w 4 if the mean photon pair number is 
small (see Appendix A). Since the transmission of 0.1 of our 
router circuitry is a technical issue, our experiment is a proof- 
of-concept demonstration of single-photon multiplexing and 
demonstrates its advantage. Based on current technology and 
theoretical analysis of the fiber coupling efficiency, it may be 
feasible to improve the transmission per router to about 95%. 
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FIG. 5: (Color online) Dependence of the photon correlation g ( ' (0) 
vs the heralding output probability P ollt put, determined in a numerical 
simulation of SPDC based HSPSs with different types of detectors 
(see Appendix D and Ref. |41 1). The heralding output probability 
Poutput is defined as the ratio of the twofold coincidence counts (be- 
tween trigger and output) and the laser repetition rate. The blue curve 
is for a 1-SPDC source with a standard bucket photon detector (D) 
with 60% detection efficiency for the trigger photon. The red curve 
is assuming a photon-number-resolving detector (PNRD) with 95% 
efficiency. Remarkably, as shown by the green curve, a 4-SPDC 
multiplexed system with standard detectors is able to outperform a 
1-SPDC heralding source realized with a highly advanced PRND. 



This number is obtained from experimentally achievable Fres- 
nel losses on each optical surface and theoretically predicted 
fiber coupling efficiency [36]. Then, the total gain of a 4- 
SPDC source with routers compared to a 1-SPDC source with- 
out router would be approximately GT w 4 • 0.95 2 = 3.61. 

Practically, the repetition rate of multiplexing is only lim- 
ited by the repetition rate of the pulsed laser (in our case 80 
MHz, but femtosecond lasers with 1 GHz are commercially 
available [37|), the jitter of the detector (500 ps in our case, 
but those with 50 ps are commercially available ll38l ) and the 
switching time of the modulators (5.6 ns in our case, but less 
than 100 ps has been shown [26 1). In addition, by employing 
state-of-the-art photonic structures, low insertion loss of the 
modulators can be achieved 11261 . It is thus possible to realize 



a gigahertz on-demand single-photon source with this scheme. 

Since in our experiment the photons are generated by 
the SPDC process pumped with a femtosecond-pulsed laser, 
which has been successfully used in multiphoton experiments 
in the past, a single photon source based on our scheme is 
suitable for high-quality multiphoton interference [31], which 
is at the heart of the quantum repeater [39| and linear optical 
quantum computation |6|. It is very interesting to note that a 
4-SPDC source such as demonstrated here based on standard 
bucket single-photon detectors, would potentially outperform 
a heralded photon source based on a highly efficient photon- 
number-resolving detector (PNRD) such as a superconducting 
transition-edge sensor (TES) |40| in specific multi-photon ex- 
periments Pffl (see Fig. [5] and Appendix D for details). This 
shows, that a 4-SPDC source would have important practi- 
cal implications for advancing existing linear optical quan- 
tum computing experiments. It is conceivable that a scheme 
which combines the multiplexing technique with supercon- 
ducting detectors as trigger detectors will be advantageous. 
Moreover, the scheme we have presented here is universal for 
any photon-pair source, including SPDC, spontaneous four- 
wave mixing, and so on. 

For the future realization of a multiplexed system with 
many levels, stabilizing the interferometers will be challeng- 
ing for bulk optics but less for integrated optics. Moreover, a 
60% transmission throughput of an integrated photonic quan- 
tum circuit has already been reported 114211431 . which includes 
the fiber-coupling efficiency and is well above the loss thresh- 
old of 1/2 mentioned above. Therefore, together with state - 
of-art micro-optics technology, (for instance, integrated pho- 
tonic quantum circuits on a silicon chip 14211431 . faster mod- 
ulators [44 1, and on-chip single-photon detectors [45 1), it is 
possible to develop a compact, scalable and nearly on-demand 
single-photon source for photon-based quantum science and 
technology following our scheme. 
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Appendix 

A. The model of an array of multiplexed heralded single 
photon sources 

In an array, m sequentially pumped HSPSs, based on 
SPDC, are coupled with m — 1 photon routers. Since the array 
of SPDC sources is pumped sequentially, the source with the 
shortest path provides the first chance for getting a single pho- 
ton, the next longer provides the next chance, and so on. The 
routers are operated in a way such that once they are activated 
by the trigger detector fired in a particular path, they block the 
outputs of all the other sources. Therefore, the source with the 
shortest path has "priority" in the array. The overall probabil- 
ity of obtaining one and only one output single photon out of 
the whole array is as follows: 



Q 1 = P 1 + P Pi 



1 — P m 

■FT' 1 Pi = Pi °— 

1 - Po 



(3) 



where Po is the probability of generating output photons 
from an individual SPDC source. Analogously, the probabil- 
ity of obtaining more than one output photon is: 



Q 



>i 



i pm 

P ° 



1-Po 



(i-Po-iV ° 



l-Pn 



(4) 



Therefore, the signal to noise ratio of obtaining one output 
single photon from this array equals that of a single SPDC 
source: 



Qi 



Pi 



= Pi 
Q>i (i-Po-Pi) P>i 



(5) 



The probability of obtaining one output single photon is en- 
hanced by a gain factor of: 



G 



Qi 

Pi 



i - PS' 

l -Pa 



1 



1-Po 



(6) 



where the approximation at the end holds for P™ <C 1. 

In any experimental setup, photon losses in the routers de- 
crease the performance of the multiplexing scheme. We now 
analyze the influence of these losses quantitatively. We con- 
sider the case where the mean photon pair number N is much 
smaller than 1 . Then, P is close to 1 and thus the gain be- 
comes G w m [as long as m(l — pj) <C 1]. Due to the cas- 
caded setup, however, using m SPDC sources means that the 
output photon has to travel through n = log 2 m path routers 
(assuming log 2 m is an integer number) and hence has an in- 
creasing chance of being lost. The output photon's transmis- 
sion through all n path routers be denoted as T. In order 
to have an advantage over a single photon source, the gain 
G has to overcompensate the loss such that the total gain, 
G to t = G T is larger than 1 . This leads to the generic con- 
dition T > 1/m. Then, G to t is increasing with m and there 
is no upper bound. In case the transmission through all path 
routers is the same, and each single router transmission is de- 
noted as Tq, the output photon's total transmission through the 



routers is T = T, 



l.i,- 



The total gain 



Gtr 



GT&mT, 



log 2 m 







(7) 
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has to be larger than 1, which leads to the condition that the 
transmission per path router must fulfill To > 1/2- Therefore, 
multiplexing is advantageous over a single-photon source if 
and only if the loss of each individual path router is smaller 
than 1/2 (under the assumption of small mean photon pair 
number). Losses in the routers do build up but are then over- 
compensated by the gain of the multiplexing. 

In our hybrid experimental scheme, two different kinds of 
routers are needed, because we employ the polarization and 
the path degrees of freedom. In an m-SPDC arrangement uti- 
lizing polarization and path, there are only ™ — 1 path routers 
(in our case, only one MZI because in = 4), and it is suffi- 
cient to have one polarization router (in our scheme, EOM3) 
at the end of the array. Therefore, every output photon has 
to travel through log 2 ? path routers (assuming log 2 ^ is an 
integer number) and one polarization router. When we denote 
the transmission of each path router by To and the transmis- 
sion of the polarization router by T po \, the total transmission 

of the hybrid multiplexed system is Th y b r = T Q 2 2 T po \. The 
total gain therefore becomes 



gT'^gt, 



hybr 



mT, 



log 2 



"T r 



pol • 



(8) 



Demanding the total gain has to be larger than 1, leads to the 



condition Tq > 



In the case of a 4-SPDC multiplexed 



2 V 2T P" 

system with high-transmission polarization routing, T po \ w 1, 
and the requirement for the transmission of the path router is 
To>i 



B. The pulsed four-photon spontaneous parametric down 
conversion sources 

We used ultraviolet (UV) femtosecond pulses of an up- 
converted mode-locked Ti: sapphire laser to pump a 2 mm 
type-II BBO crystal at a wavelength of 404 nm with a pump 
power of 1400 mW, a pulse duration of 180 fs, and a repe- 
tition rate of 80 MHz. The UV pump power was adjusted 
by varying the up-conversion efficiency via moving the up- 
conversion crystal in or out of the focus of the lens for the fun- 
damental pulse. The UV pulses successively passed through 
two BBO crystals to generate two polarization entangled pho- 
ton pairs (photons 1 and 2, and photons 3 and 4) in spatial 
modes a and b, and c and d via type-II spontaneous paramet- 
ric down-conversion respectively, as shown in Fig. 2(a). Af- 
ter proper filtering, maximally polarization entangled states of 
the form ^(|#> a(c) 



\V) 



b(d) 



l^)a(c) I #)b(d)) were gener- 
ated. 9 is the phase difference between horizontal and vertical 
polarization due to birefringence in the crystal, which could 
be corrected by placing a half wave plate (A/2) oriented at 
45° and 1-mm compensation BBO crystals in each photon's 
path [23 1 . By tilting the compensation BBO crystals, we set 
the phases to be for both sources, and hence we obtained the 
two polarization-entangled states described in the text. 

To exactly define the spatial and spectral properties of the 
emitted four photons, we coupled each photon into single- 
mode fibers and filter them with 3 nm (full width at half max- 
imum, FWHM) bandpass interference filters (IF) for photons 



1 and 4, and 1 nm IFs for photons 2 and 3. In the HOM in- 
terference experiment, we used 3 nm IFs for both photons. 
The photons were detected by single-photon counting mod- 
ules (SPCMs), and the coincidences were registered with a co- 
incidence logic. At the pump power of 1000 mW, the average 
trigger rate is about 500 kHz and the mean photon pair num- 
ber generated per pulse is approximately 0.062. Note that the 
coupling efficiency of the sources from BBO 1 is about 13%, 
and that from BBO 2 is about 7%. This difference eventually 
gives the difference in the correlation functions and the count 
rates of each source. For clarity, the 1-SPDC data shown in 
Fig. H were obtained by averaging the g 2 (0) functions and the 
counts of the four individual 1-SPDC sources. The 2-SPDC 
(4-SPDC) data were obtained by averaging the g 2 (Q) func- 
tions and summing the counts of two (four) 1-SPDC sources. 



C. The ultrafast photon router 



We used a MZI to realize the ultrafast photon router. The 
performance of this router strongly depends on the phase sta- 
bility of the MZI. To achieve that, we built the interferometer 
in an enclosed box made by acoustic isolation materials in or- 
der to stabilize the phase passively. Additionally, an active 
phase stabilization system has also been implemented. This 
was accomplished by using an auxiliary beam from a power- 
stabilized He-Ne laser counter-propagating through the whole 
MZI with a little transversal displacement from the signal 
beam, which picked up the phase fluctuation of the interfer- 
ometer. After passing through the interferometer, the inten- 
sity fluctuation of the He-Ne laser beam was measured with a 
silicon photon detector and the signal was fed into an ana- 
log proportional-integral-derivative (PID) regulator. A ring 
piezotransducer attached to one of the mirrors in the MZI was 
controlled by this PID regulator and compensated the phase 
fluctuation actively. The insertion loss of this router and fiber 
delays was about 90%, which was mainly due to the single- 
mode fiber coupling, lossy narrow-band interference filters (in 
spatial modes b" and c", 1 nm FWHM, about 70% peak trans- 
mission) and Fresnel loss on various optical elements. The 
most straight-forward step to improving the transmission of 
the router is to use high peak transmission and broadband in- 
terference filters. With a 95% peak transmission and 3-nm 
FWHM interference filter, we have measured the total trans- 
mission throughput of the router, including the delay fibers 
and correlation function measurement setup, to be more than 
25%. Note that one can engineer the interactions between the 
pump and nonlinear crystals and match the group velocities of 
the pump and down-converted photons [ 46-48 1 . By employ- 
ing this novel technique, one can generate pure single pho- 
tons without using narrow-band filters, which potentially al- 
lows us to further enhance the transmission throughput to be 
more than 50%. 
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D. Comparing 4-SPDC to a 1-SPDC with different types of 
photon detectors 

In order to estimate the required performance of HSPS in 
practical applications, we implemented a numerical simula- 
tion of the quantum optics circuits. Thereby, every photon 
mode is represented in a Fock space from to 4 excitations, 
and all unitary operations are connected Hamiltonians ex- 
pressed within the full Hilbert space. Hence, all higher order 
photon terms are intrinsically taken into account, up to the nat- 
ural cutoff from the Fock space. Through numerical simula- 
tion we studied the output probability for the HSPS versus the 
observed second-order correlation function for zero time de- 
lay, <? 2 (0) = ?J a \2 , for various pumping strengths e 2 of the 



SPDC operator, -ffsPDC = e(a^ +ba), where a and b are the 
two photon modes. These results are shown in Fig. [3] Inter- 
estingly, the 4-SPDC source based on 60% efficiency bucket 
detectors (for detecting trigger photons) and optical transmis- 
sion throughput of 95% per router, which should be reason- 
able to reach with an optimized optical system, can outper- 
form a HSPS built from 1-SPDC heralded with the currently 
best possible detector technology, such as a PNRD system 
(for detecting trigger photons) with total detector efficiency 
of 95%. Note that we assume the coupling efficiency from the 
source to the detectors to be 70% for both cases, because that 
is independent of the detectors. 



